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Using the modified internal shock wave model we fit the gamma ray burst (GRB) light 
and spectral curves of 30 GRBs observed with BATSE. From the best fitting we obtain 
basic parameters of the relativistic shells which are in good agreement with predictions 
given earlier. We compare measured GRB parameters with those obtained from the 
model and discuss connections between them in the frame of the physical processes 
laying behind GRB events. 
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1. Introduction 

The quest on resolving a Gamma Ray Burst (GRB) event consists of finding an 
explanation for several parts e.g. spatial distribution of the event, afterglow, spectral 
and light curve, collimation, etc. In the past decade the GRB phenomenon has 
been thoroughly investigated both observationally see e.g. PQ-P3] and theoretically 
see e.g. [15]-[23]. However, a mystery of GRB phenomena lies in its heart, where 
the central engine ejects material with relativistic energies and velocities. Due to 
the high optical depth of the expanding material in the first phase of the GRB 
event, structural observations of the central engine, that is located near the core 
of the progenitor, are not possible. The only information that one can obtain from 
the observations, in the first minute of a GRB event, is the temporal variability 
of the 7-ray light curve. This usually shows strong, but short fluctuations of the 
energy output with a typical time-scale of the order of milliseconds to a couple of 
seconds [24]. The numerical simulations of [25] and [26] for the early phase of the 
explosion, revealed that 7-ray light curve pulses replicate temporal activity of the 
'inner engine'. This can give information about the connection between observations 
and physical processes occurring in the GRB core. 
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In order to make conclusion about the physical processes in the GRB core one 
should assume the mechanism of GRB origin. The mechanism of pulse creation in a 
GRB light curve is proposed to be connected with mutual interactions (collisions) 
between faster and slower spreading shells see e.g. [22], in the so called internal shock 
model. Here we accepted slightly modified version of the model with the additional 
assumption of the non-zero density environment and also with different treatment 
of the slower material (shell) accumulated at some distance from the GRB engine. 
In a previous paper [27] , where we considered interaction of an incoming faster shell 
with dense barrier, we demonstrated that the model is able to reproduce and fit 
well the observed GRB light curves. 

In this paper we apply the model to a sample of 30 GRB light curve pulses 
observed with Burst and Transient Source Experiment (BATSE) in order to discuss 
the physical parameters of relativistic shells. 

The paper is organized as follow: in §2 we give a brief description of the model 
and considered assumptions; in §3 we describe the observational sample and fitting 
procedure; in §4 we discuss our results and finally in §5 we outline our conclusions. 

2. The model and assumptions 

Here we describe the physical scenario and some approximations used in the model 
of GRB light and spectral curves. First, we assume that the GRB engine ejects 
an amount of relativistically expanding material, that spreads isotropically from 
the center of explosion. The material is subsequently ejected from time to time 
depending on the central engine activity. Here we will not consider the nature of 
GRB progenitors, i.e. for the model it is not relevant whether a GRB event is 
originated in the collapse of a massive star [28] or in the process of merging of two 
compact objects - neutron stars or black holes [29], [30]. The most important for the 
model is the assumption that the ejected material can mutually interact or interact 
with the surrounding environment. The ejected material is probably irregular in 
nature, with different initial parameters (mass M e j and Lorentz factor Fo). 

During the expansion of the ejected material (closer to the GRB core), a slow 
moving material is followed by a fast moving one, thus the faster moving material 
will overtake the slower material and plunge into it. This interaction produces the 
relativistic shells and shock waves that can accelerate particles to very high speeds. 
The new formed relativistic shells are probably with different velocities and could 
collide mutually producing the observed GRB light curve pulses. The observed 
GRB light curves show vast variety of pulses, ranging from very intense ones to 
those almost equal to noise, and from symmetrical to highly asymmetrical ones. 
To model different observational light curve pulse profiles, we consider here the 
modified standard internal model, assuming that the accumulated matter is able 
to form a slow moving barrier. This, allowed us to generate more diversity in light 
curve pulse profiles. 

In the shell interaction processes, the mass of a newly created shell is approx- 
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imately the sum of the two colliding shells and its Lorentz factor is smaller than 
it was in the faster shell. It could be additionally reduced by surrounding medium 
or further shell interactions, thus producing the accumulation effect. Depending on 
the density and width of such created dense barrier, further collision with it pro- 
duces wider/thinner or higher/lower pulses. Also, considering the energy density of 
expanding shells and barriers, one can get more symmetrically shaped light curve 
pulses in case of the shells with high energy and low energy barriers, and high asym- 
metrical pulse shapes in the opposite case. This is the main mechanism that we use 
to explain high temporal variability of the observed 7-ray light curves. From the 
scenario described above, it follows that the selection of shell parameters can have 
a random distribution in a given interval of values. 

Similar as in [25] , we consider that the ejected GRB material is organized in 
an ultra-relativistic flow of well defined and collimated shells with random initial 
energy. In contrary to [22] and [25], where null density hypothesis was used, we 
assume here that the surroundings regions around the GRB central engine consist 
of at least small number of particles with densities uq > 1 cm~ 3 . This allows us 
to analyze the density of the moving shell and to model the density distribution 
and shape of the barrier in a specific way (here we use Gaussian function, see Eq. [5] 
further in the text). The assumption of no 7^ around the central engine seems to 
be generally valid. If there are some kind of ejections from the central engine, one 
can expect an amount of scattered material distributed to the surrounding region, 
as e.g. in the collapsar model where a central star is Wolf-Rayet type. Such star 
ejects a huge amount of material during its final stage, therefore one can expect non 
zero density in its environment. 

To describe the evolution of a relativistic shell, we adopt a phenomenological 
model based on the [TH] that presents a system of the first order of differential 
equations (where the distance R, Lorentz factor T and mass m s of the shell are 
included, see [27]): 

dR 



^ = cVf^i[r + Vf^i\, (i) 



dm s M ei + 2{l-€)Tm s +Zm s , w 



dm s , „,R 2 f dR dY\ 

27rnm p (l - cos 9) — [3T— - 2R— , (3) 



dt pv 1 T 3 V dt dt J ' 

where the parameter £ takes values from in case of adiabatic expansion to 1 in 
fully radiative case. M e j and 9 are the initial mass and collimation angle of the shell, 
n is the number density and m p is the proton mass. Eqs. (|T|) - (J3j) are derived for 
an observer reference frame, and they have to be solved simultaneously, together 
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with the density equation. The initial values of parameters and variables are highly 
dependent on physical properties of the shocks. 

The density of the barrier cerated from accumulating decelerated shells (emitted 
by the central engine) could be described with (see [3Tj): 



where n% and n\ are number densities after and in front of the shock, K2 is the ratio 
of the specific heat for the shocked fluid and 72 is the Lorentz factor of the shocked 
fluid. 

This equation gives a connection between density of perturbed and unperturbed 
material. In the case of the ultra-relativistic expansion, the ratio of the specific heats 
has a constant value of n = 4/3, then Eq. [4] can be reduced to 2a. = 472 + 3. Also, 
in the relativistic regime, the Lorentz factor of the shell F is directly proportional 
to the Lorentz factor of the shocked particles 72 [3T| . 

In case of the collision of relativistic shells, a slower interacting shell (or barrier 
in our model) presents a density perturbation in surrounding media for the incoming 
faster shell. The barrier (or accumulated material) at a distance R c from the central 
engine should have a distribution of the density along the path of penetration of 
the faster shell. This can be included in calculation, and here we assume a Gaussian 
density distribution as: 



where a and b describe the Gaussian intensity and width of the barrier and hq is 
the density of the surrounding region. i?o is the initial position of the faster shell, 
F is the Lorentz factor of the shell as used in Eqs. (1-3). Note here that in Eq. 
[5] (as it was mentioned above) uq > 0, and it is a crucial difference between our 
and the standard internal shock (IS) model. In the standard IS model, the density 
of the interacting shells is fixed by the mass loss rate from the central engine, the 
Lorentz factor of the shells and the distance from the source where the collisions 
take place. It is not directly related to the density of the external medium. As a 
result, even in the absence of an external medium, prompt emission will result from 
shocks taking place in the material ejected from the source. Indeed there is some 
evidence of bursts occurring in very low density environments which have a prompt 
emission but no detectable afterglow, but as we mentioned above, one could expect 
that no > (especially close to the accumulated material) in the central engine 
surroundings. 

Generally, one can expect that the ejected material in form of the shell can collide 
with the ISM which has a certain density distribution. In such a highly relativistic 
physical system the relative motion of charged particles of the ISM can generate 
an intense magnetic field in the reference frame of the moving fluid. We calculate 
the magnetic field in a similar way as in [32 , by assuming that the energy of the 




n 2 _ K2J2 + 1 
m K2 - 1 



(4) 
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magnetic field is a certain fraction, of the total energy of the relativistic shell. 
In the comoving reference frame the magnetic field is calculated as: 



B = y87r£ fc n rm pC 2 (4r + 3) \?f^j (l + a ■ eA^)' J , (6) 

where the variables used in this equation are same as in the Eqs. [l][5j 

The emission mechanism of shock waves is mainly based on the synchrotron 
radiation, but for higher energy bands additional flux may be gained by the Inverse 
Compton (IC) radiation [22]. In the first approximation we neglect the IC radiation. 

We calculate the intensity of the radiation emitted by particles in the relativistic 
shell using the formulae given by [SHI ■ Then the total emitted flux can be calculated 
as e.g. in [33]. Note here that an expanding shell contains relativistic electrons and 
baryons which contribute to the synchrotron radiation. However, taking into account 
the difference in velocities of these constituents, one can neglect the contribution of 
baryons to the total emitted flux. Then in the comoving reference frame the total 
flux is given as: 

K = A- l^ +1) F{v / V ' syn )d le (7) 

where A, j emin and j ema x are: 

_ V3e 3 B' m s _ _ (p - 2) mp 

m e c z m p (p — 1) me 

and 

/•OO 

F(x) = / K 5/3 (x) dx, 

J X 

where, K 5 / 3 is the Bessel function of the second order and v syn = 3j^eB /47rm e c 
is the critical frequency of the synchrotron radiation. 

In Eq.[7]we neglect the effects of the surface curvature of an emitting shell, since 
it has a small influence on the pulse shape. 

In the case of ultra relativistic shells the cooling time is much shorter then 
the dynamical (time of expansion) (see [35]). This is particulary interesting in the 
gamma phase of explosion where shells are interacting with each other. 

With the model described above we are able to simulate collisions of relativistic 
shells in the first phase of a GRB event, which produce the peaks in the light curve 
[27] . This model has been used to fit the light and energy curves of the sample of 
30 GRBs taken from the BATSE database. 
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3. Model vs. observations 

In [27] we demonstrated that the model is able to reproduce (simulate) the observed 
light curves of GRBs. Moreover, the model can properly fit the observed light curves. 
Here we selected 30 GRBs from the BATSE database and fit them with the model 
in order to find the basic parameters of interacting shells. In this section we describe 
the selected sample and the fitting procedure. 

3.1. The sample 

From a large BATSE database (3rd channel, E — 100 — 300 keV, for the light curve) 
we select a sample of 30 GRBs using following criteria: (i) GRBs have isolated light 
curves with the clear peak maximum. For proper application of the model we avoid 
the complex pulse profiles. @ (ii) we avoid small pulses because of their low temporal 
resolution; (iii) we include in the sample as much as possible different GRBs (long 
and short lasting, with different profiles, different intensity and different profile 
asymmetries). 

In Table[T]we give a list and basic parameters of selected GRBs. The parameters 
in the table are (from the first to the last column respectively): the Full Width at 
Half Maximum of the intensity of light curve pulse (FWHM), the time of peak 
intensity for the observed pulse (£ pea k), the total duration of the pulse from the 
beginning to the end of its lower tail (At), the maximal intensity (J m ) of the pulse 
measured in erg/cm 2 sHz, and the asymmetry indicator (w) calculated as a ratio of 
half-halfwidths before and after the maximum. 

Fig. [1] shows statistics of the measured parameters given in Table [TJ As it can 
be seen in Fig. [1] the values of observed parameters do not follow the Gaussian 
distribution. Due to the relatively small number of GRBs in the sample, here we 
are not able to discuss the power law indices of parameters (for detailed studies of 
these GRB properties see [24], [36]). 

3.2. Fitting procedure 

In order to follow the relativistic shell evolution and collision, we consider the case of 
only one shell expanding from the GRB core. It propagates through the surrounding 
media, which can contain a barrier with the mentioned Gaussian profile. We fit with 
our model light and energy curves of GRBs given in Table [T] In order to find the 
best fitting we vary the parameters of the faster shock and barrier. We consider 
the following parameters as free: the Lorentz factor To, the total initial ejected 
mass of the shell M e j, the density of surrounding media no, the opening angle of 
the jet 9 m , the distance of collision R c and the parameters a and b which describe 
the shape (height and width) of the density barrier. Additionally, we assumed that 



a The process of pulse creation is stochastic in nature and may result in a complex pulse shape, 
e.g. it is often observed that two or more pulses are superposed; 
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Fig. 1. Histograms of basic GRB characteristics for the sample of GRBs presented in Table[T] On 
the Figures (a) to (d) we place distributions of (a) the Full Width at Half Maximum of the intensity 
of light curve pulse (FWHM), (b) the time of peak intensity for the observed pulse (i poa k), (c) the 
maximal intensity (J m ) of the pulse and (d) the asymmetry indicator (w) 



the barrier can move, thus we also put as a free parameter the Lorentz factor of 
the barrier IV In Fig. [2] (left panels) we show the best fit of three isolated pulses 
with different shapes: GRB000508, GRB911104 and GRB911117. The light curves 
of these GRBs do not have a standard form, i.e. the shape of pulses does not always 
follow the FRED (Fast Rise Exponential Decay) behavior. As one can see from Fig. 
[2] the shapes of the light and energy curves can be very well fitted with the model. 
In right panels we show the best fit of the averaged spectral energy distribution 
(ASED) taken from all four BATSE channels. Although the data for measured 
counts in the energy channels are with large uncertainties, the fit of ASED can be 
used for the confirmation of the validity of the GRB light curve fit, since the same 
parameters are used to fit both curves. 

The coefficients £, £ e and ^JB, determine the radiation efficiency for expanding 
relativistic shell. The synchrotron radiation is directly proportional to intensity of 
the magnetic field, as well as the energy of the electron component of plasma. 
Consequently, if one puts high values for those three components it will increase 



b £ describes fraction of total shell energy that has been converted into radiation, £ e fraction of the 
total shell energy devoted to the electron plasma component and fraction of total shell energy 
contained in the magnetic field. 
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the intensity of radiation and also the intensity of light curve pulses. However, the 
conversion from the kinetic to radiative energy is the process with relatively low 
efficiency (see e.g. [37]), therefore the above coefficient must not exceed 10 to 20% 
[22] . They are expected to be smaller than 0.2. In order to reduce the number of 
free parameters in the fitting procedure we fixed them to £ = 0.1, £ e = 0.2, £j = 0.2 
(see also [38]). Moreover, the different values of these parameters will mainly affect 
the intensity of GRB light curve, but not the shape which is mainly considered in 
the fitting procedure. 

Also, the distribution of electrons in the shell follows a power law function where 
the index of the electron distribution p usually takes a value between 2 and 3 [39] . 
In the fitting we fix it to p = 2.5 which well reproduce the obtained values from 
the fits of the GRB light and energy curves. In the fitting procedure we used x 2 
minimization. 

4. Results and discussion 

The parameters of the shells and barriers for 30 GRBs obtained from the best fittings 
are given in Table [U The distribution of parameters are presented in Fig. [3] As it 
can be seen from Table [5] (also see Fig. [3$ the most of GRB light curve pulses in our 
sample indicate an initial Lorentz factor around 95, with the mean value of about 
93. Cases with Tq < 90 could be explained if we consider the physical mechanism 
of the collisions. In the first place, the initial value of the Lorentz factor directly 
determines the initial energy and velocity of the incoming shell. So, for those shells 
(higher Tq) produced light curve pulses will be more intense and short lasting. On 
the contrary, smaller Tq will produce low intensity and long lasting pulses. This is in 
a good agreement with the observations of the pulse width - luminosity correlation 
found in [36] . 

The total released isotropic energy of a GRB lies in the range from 8 x 10 47 
erg (GRB 980425 associated with supernova SN1998bw, see e.g. Ref [40], [41]), to 
2 x 10 54 erg (GRB 990123 Ref. 0). Note here that using the conical jet model 
|42) found that the gamma-ray energies are clustered in a narrow interval around 
5 x 10 50 erg. The mass of the ejected shell in the most GRB light curve pulses is 
in the interval of ~ 10~ n — 10 -10 M sun , with a mean value of 10 _10 M sun . These 
could be used to calculate the initial energy of a particular shell in the moment of 
the ejection, giving values in the range from 10 45 to 10 47 ergs. The sum of energies 
of all ejected shells during the GRB event can indicate the total energy released by 
the central engine in the particular GRB. 

Additionally we assume that the barrier is moving, with the Lorentz factor Th- 
in the first approximation we take that its velocity is constant until the collision. 
This parameter has multiple influence mainly on the intensity and width of the 
created light curve pulse, as well as on the FRED pulse shape. For the faster barriers 
interaction is long-lasting, but with a low intensity. On the other hand, if the barrier 
is moving significantly slower than the incoming shell, it will cause a more intense 
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Fig. 2. The light (dashed line - left panels) curves and the averaged spectral energy distribution 
in four BATSE channels I ch:(20-50)keV, II ch:(50-100)keV, III ch:(100-300)keV, IV ch:>300keV 
(full circles - right panels) for GRB911104, GRB911117 and GRB000508 (from top to bottom, 
respectively) fitted with the model (solid lines - left panels and dashed lines right panels). 



interaction followed by the short-lasting and very strong pulses. The FRED shape 
is more dominant in the former case. 

The density of the ISM in the region around the central engine is assumed to be 
homogeneous (s = 0), with a density approximately an order of magnitude higher 
than at the distances where afterglow starts (see for example [43]). We obtained 
values from several to few tens of particles per cm -3 (see Tabl^2|) that is in the 
expected range n ~ (10° — 10 3 ) cm~ 3 . The distribution of n (Fig. [3jl) for our 
sample has a bell-like shape with the maximal value around 60 cm -3 . This values 
is more appropriate to hypernova then to a merger scenario, because in the merger 
scenario one can expect significantly smaller densities |36) . 

Similar distribution is obtained for opening angle of relativistic shells 9 m , which 
is considered to be constant during the evolution. The most probable values are 
around 0.05 radians (~ 3°) (see [42] and reference therein). Existence of this par- 
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Fig. 3. Histograms for obtained parameters from the best fit: To, M e j, Tb, n , 9 m , R c and AR 
(a-g panels, respectively). 

ticular value is determined by physical processes in the vicinity of the GRB central 
engine. Namely, if one take a higher value of 8 m , the resulting pulse is broader and 
a slow decay feature is much more visible than in the case of smaller angles, where 
the pulse is thinner and has a symmetrical shape. That is in a good agrement with 
engaged physical processes during the shell expansion. 

We suppose that the shell interaction occurs mainly close to the GRB engine, 
(distance R c ~ 10 14 cm), as it was proposed in the internal shock scenario (see 
[22]). This parameter has very small influence on the shape of GRB pulses, but has 
an influence on the intensity of pulses. Pulses produced in a collision closer to the 
GRB engine are more intense than ones originating at larger distances. 

The parameters of the Gaussian of the density distribution which describe the 
structure of the barrier, the width at the half maximum b and the density in the 
central part a, both influence the shape of GRB light curve pulses. They are trans- 
lated into appropriate variables, number density rib and width AR of the barrier, 
respectively, using the connections given by Eq. ® and [TUJ 



n b = n (l + a) 



(9) 



AR = 2b\ fin 

V a — 1 



(10) 
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Note here that in the case of dense barrier one can expect n^jriQ >> 1 (see 
also Table [2]), consequently a » 1 and Eq. [10] can be rewritten as AR = 26V In 2. 
The influence of the barrier parameters on pulse profiles is following: a barrier with 
narrow width and high number density will have a strong interaction with fast shell, 
producing symmetrical and intense light curve pulses. On the contrary, a barrier 
with larger width and lower number density will cause small intensity pulses with 
high asymmetries. 

In general, comparing the obtained values of parameters (Table [2]) for different 
GRBs, one can conclude that there is no significant difference between them even 
when the shapes and durations of GRB pulses are different. This suggests that the 
nature of GRBs is similar and that there should be no big difference between the 
physical conditions of GRB progenitors (see [3]). On the other hand, the barrier 
density distribution can differ from the Gaussian, that is assumed here, and it may 
reflect the values of basic parameters. But in any case, one can expect that the 
density distribution of the barrier has to be taken into account in the shock model. 

4.1. Connection between the shell parameters 

In order to find physical meaning of the obtained parameters, we explore correla- 
tions between them. In Table [2] the parameters are divided in two groups: one that 
describes a shell and other connected with the barrier. We expect that parameters 
from those two groups are not in correlation, since they are independent. The re- 
sults are presented in Figs. |4] and 03 where we separately denoted the long (t>2s) 
and short (t<2s) GRBs with open and full circles, respectively. As one can see in 
Figs. 0] and [5] there are no strong correlations between these quantities, but only 
in some cases, there is a slight connections between different parameters, as e.g. 
6 m vs. Tq (FigHk), M c j vs. Tq (Fig. HJd). In the case of the barrier parameters a 
slight correlation can be found in AR vs. nf, (Fig. [5k,). Additionally, we examine the 
correlation which may be established between Lorentz factors of the incoming fast 
shock wave Tq and the moving barrier Tf,, presented in Figure EJ). 



a : iof b - 

- i 8 j - 

'„ ' "". s 2 ° . ; • 

° ' ■ d p i . " ° • . > . -J 

70 80 90 100 110 120 130 60 70 80 90 100 110 120 130 
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Fig. 4. Parameter dependance: 8 m vs. To (a) and M e j vs. To (b). The long (t> 2s) GRBs are 
denoted with full circle and short with open circle. Dashed line presents the border above which 
there is no parameter values for this sample of GRBs. 
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Weak correlations or non-correlations can be noticed in each case particulary. 
For example, with larger opening angle of the shock wave 9 m , the volume of the 
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shell increases, that causes smaller initial Lorentz factors of the shell and vice versa 
(see Fig. [4^,). For the certain shell energy, the increase of the space angle 8 m causes 
the increase of the ejected mass M e j, and that causes the decrease of the initial Tq. 
This is in a agrement with the trends presented in Figured, where one can see that 
for a higher value of M e j, the initial Lorentz factor tends to be smaller. Of course, 
one can not expect to see obvious confirmation for this conclusion with such small 
sample of examined GRBs, but rather signs of trends for mentioned dependence. 
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Fig. 5. Parameter dependance: AR vs. ni, (a) and To vs. Ff, (b), with r designating the correlation 
coefficient. The notation for short and long pulses is same as in Fig. 4. 



In the case of the barrier parameters, there is an indication of the connection 
between the particle number density rib and width of barrier AR as it is shown in 
Fig. [5^, as AR ~ l/n&. Physically one can expect such situation, i.e. for a broader 
barrier, the density of barrier tends to be smaller and vice versa. 

In Fig.Eb the r vs. r& is presented. There is no global correlation between these 
Lorentz factors, but taking into account only short pulses there is some indication 
that for a faster barrier the Tq of the shell is higher. The dashed and solid line in 
Fig. EJd present the linear dependance of the given parameters, indicating that the 
incoming shell must have higher Tq than the slower moving barrier, providing just 
enough necessary conditions for event of collision to happen. 

4.2. Connection between shell parameters and observed pulse 
parameters 

Additionally, we explore possible correlations between parameters obtained from 
fitting the light curves and measured ones (given in §3.1). In Fig.[6]we plot measured 
values against Tq ■ M e - } , which is proportional to the energy of the incoming shell. 

It is interesting that the pulse intensity for small energies (Tq ■ M c - S <0.2) shows 
nearly linear trend with the energy (see Fig. [6fc), but for Tq ■ M e j > 0.2 this trend 
is not present. 

In Figs[7K and b we present At vs. Tq, Tj,, respectively and in Fig.[7t FWHM vs 
AR. It is obvious that for a faster barrier and shell the interaction will be shorter. 
Also, for a broader barrier we obtain long lasting (wider) pulses as it is shown in 
Fig. [7b. There is a correlation between FWHM and AR, with correlation coefficients 
r = 0.61 for short and r = 0.84 for long lasting GRBs. A linear relationship between 
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Fig. 6. Parameter dependance: FWHM, At, J m and w vs. To ■ M e \ (scaled to the Maxpo ■ M e j]), 
panels a - d, respectively. The notation for short and long pulses is same as in Fig. 4 and 5. 



FWHM and AR is present as FWHM=a-Ai? + b (where a = 0.54 and b = 0.22 for 
short and a = 1.2 and b = 0.024 for long lasting GRBs, see Fig. [Th). 
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Fig. 7. Correlation of pulse and shell parameters. Case of Af vs. To, T;, (panels a, b), and FWHM 
vs. AR, (panel c). The notation for short and long pulses is same as in Fig. 4 and 5. 



We found some connections between n& vs. AT and n& vs. J m , as presented in 
Fig. [8] This is expected, since for a barrier with higher density there are more intense 
pulses (Fig. EK). The smaller values of density produce broader light curve pulses 
(Fig. [8b). It can be seen from Fig.EK that short pulses have mainly smaller densities. 
This can be easily understood analyzing the physical processes in the moment of 
collision. When the particle number density is higher more particles are taking part 
in the interaction, producing a higher gamma-ray flux. But if the decreases the 
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interaction is prolonged, caused by the stretching ol the barrier material, so long 
lasting pulses are produced. 




20 30 40 50 60 
/„, x 10~ 27 [erg/cm 2 s Hz] 



Fig. 8. Correlation of pulse and shell parameters. Case of n b on AT (a) and n b on J m (b). 



5. Conclusions 

In this paper we extracted the basic parameters of internal shock waves during the 
first phase of a GRB event by fitting 30 observed GRBs (from BATSE database). 
To fit the observed GRB light curves we used the modified internal shock model 
given in |27j , assuming the collision of a fast shell with a slow moving (with respect 
to the velocity of the shell) barrier. We analyzed the obtained parameters in order 
to find physical processes behind the GRB origin, and we came to the following 
conclusions: 

(i) Relativistic shell parameters obtained from the fitting of GRB light curves 
are in a good agreement with expected ones and also with estimations given earlier 
by other authors mentioned throughout the text. 

(ii) The obtained values of internal shell physical parameters for GRBs with 
different light curves are in the short interval, showing that the physical processes 
behind the GRB creation are similar, i.e. there should be the ejected mass that 
collides with surrounding regions - or accumulated slow moving material. 

Also, we analyzed possible connections between parameters obtained from the 
best fitting of GRB light curves with measured ones. From this analysis we can 
conclude: 

(i) There is no correlation between parameters obtained from the best fitting 
(Figs 4 and 5), only some indication that long GRBs have higher values of Lorentz 
factor, and we found a slight trend between Lorentz factor of the shell and moving 
barrier for short pulses. 

(ii) There is a correlation between the intensity of pulses and the energy density 
of the shell only for a low energy pulses (r • M e j < 0.2, see Fig. 6c). 

(iii) The FWHM of GRB light curve pulses is in the correlation with the width 
of the barrier. Using this we give a relation between FWHM (that can be measured 
from observed light curves) and AR that is a parameter of the model (see Fig. 7b). 

Finally, we can conclude that the modified internal shock model that assumes 
also the barrier density distribution can well describe the first phase of the GRB 
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origin. Moreover, the obtained parameters for the internal shocks and barriers well 
fit the physics of the GRBs. 
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Table 1. Parameters of the selected GRB light curves. From the first 
to the last column the following parameters are given: Full width at 
half maximum (FWHM) , time of peak intensity t pea k , duration of the 
pulse At, pulse intensity (J m [crg/cm 2 sHz]), asymmetry indicator w s . 



GRB 


FWHM [sj 


At pea k[s] 


At[s] 


7 r.,in — 27i 

J m [xl0 




GRB910629 


0.55 


0.4 


1 


13 


0.6 


GRB911104 


0.85 


0.5 


1.3 


13 


0.5 


GRB920715 


0.5 


0.45 


1 


3.2 


0.4 


GRB920720 


0.75 


0.6 


1.4 


40 


0.7 


GRB920808 


0.5 


0.25 


0.7 


11.3 


0.4 


GRB920811 


0.2 


0.25 


0.3 


3.5 


0.7 


GRB920830 


2.75 


1.9 


6 


7.8 


0.6 


CJRRQ90Q1 9 


0.4 


0.58 


1.2 


5.4 


0.7 


GRB920924 


0.4 


0.4 


0.8 


4. 


0.3 


GRB921021 


1.8 


1 


4 


5.1 


0.4 


GRB921207 


1.4 


1.4 


3 


60. 


0.6 


GRB921208 


1 


0.85 


2 


2.25 


0.8 


GRB921222 


1 


0.85 


2 


2.25 


0.8 


GRB950129B 


0.6 


0.65 


1.22 


3.7 


0.7 


GRB950211B 


1.5 


0.95 


3.5 


6.5 


0.4 


GRB960111 


1.2 


0.5 


2.6 


9 


0.5 


GRB960207 


0.33 


0.25 


0.6 


9.2 


0.8 


GRB960229 


0.55 


0.35 


1 


5.4 


0.8 


GRB960311 


0.4 


0.4 


0.7 


6 


0.6 


GRB960409 


2.4 


2 


5 


10.5 


0.7 


GRB960418 


0.9 


0.8 


1.5 


4.7 


0.4 


GRB960524 


0.55 


0.4 


1.2 


6.2 


0.8 


GRB960528 


2.2 


0.55 


4.2 


4.2 


0.6 


GRB960530 


5.5 


3 


11 


5.7 


0.6 


GRB960613 


1.6 


1.8 


4 


6 


0.3 


GRB960617B 


0.7 


0.55 


1.1 


3.3 


0.8 


GRB970424 


0.4 


0.3 


0.7 


3.3 


0.5 


GRB991105 


0.65 


0.45 


1.6 


6 


0.4 


GRB991213 


1.4 


1 


2.4 


4 


0.4 


GRB000107 


0.5 


0.45 


0.7 


4.7 


1. 



January 24, 2012 1:12 WSPC/INSTRUCTION FILE ssimic5 



18 S. Simic and L.C. Popovic 



Tabic 2. The internal shell and barrier parameters obtained from the best fitting of light and energy 
curves for the sample of 30 GRBs. In the last two rows we put the mean and appropriate standard 
deviation of parameters. 



parameter 


r 






n 


9m 


Rc 


AR 


n b 


units 


- 


■10- n [M sun ] 


- 


[cm d ] 


[rad] 


■10 14 [cm] 


■10 13 [cm] 


•10 4 [cm -3 ] 


ORRQ1 0fi9Q 


109 


5.5 


50 


43 


0.05 


2.0 


5.1 


43. 




111 


10 


43 


63 


0.06 


2.5 


5.1 


56.7 


ORR92071 5 


103 


1.3 


51 


33 


0.06 


2.2 


5.8 


33. 


CRR920720 


90 


25. 


49 


70 


0.06 


1.5 


5.8 


105. 


CJRR920808 


110 


8. 


50 


33 


0.07 


1.9 


7.8 


8.3 


ORRCWORI 1 

VjllDjilUOl J. 


111 


1.1 


65 


65 


0.06 


1.2 


2.2 


58.5 


ORR920830 


71 


13.5 


45 


37 


0.1 


3.5 


25. 


1.1 




100 


10. 


74 


70 


0.05 


2.0 


3.5 


4.2 


GRB920924 


125 


2. 


50 


50 


0.04 


2.5 


3.8 


150. 


GRB921021 


95 


5. 


65 


20 


0.08 


3.5 


21.7 


0.6 


ORRQ91 907 


73 


53. 


63 


75 


0.1 


2.5 


10. 


150. 


GRB921208 


71 


3.8 


49 


70 


0.07 


1.9 


7.3 


4.9 


GRB921222 


87 


2.5 


70 


20 


0.06 


2.1 


7.8 


6. 


GRB950129B 


84 


4. 


70 


20 


0.05 


1.8 


7.5 


100. 


GRB950211B 


73 


18. 


56 


90 


0.05 


1.7 


8.3 


4.5 


GRB960111 


95 


10. 


68 


45 


0.08 


1.9 


7.7 


1.8 


GRB960207 


109 


6.8 


63 


57 


0.06 


1.5 


3.2 


39.9 


GRB960229 


98 


5.2 


61 


75 


0.057 


1.7 


4.8 


9.8 


GRB960311 


110 


13. 


77 


19 


0.04 


1.9 


5.2 


9.5 


GRB960409 


75 


18. 


52 


110 


0.1 


5.5 


30. 


5.5 


GRB960418 


75 


3. 


48 


50 


0.09 


2.7 


10.7 


50. 


GRB960524 


87 


2.7 


65 


50 


0.09 


2.1 


8.8 


335. 


GRB960528 


91 


5.5 


56 


43 


0.06 


2.0 


13.3 


2.2 


GRB960530 


62 


40. 


46 


31 


0.055 


3.8 


30.1 


0.6 


GRB960613 


87 


13. 


63 


65 


0.04 


2.4 


10.2 


3.9 


GRB960617B 


103 


3. 


87 


50 


0.07 


2.3 


8.3 


150. 


GRB970424 


99 


0.9 


53 


35 


0.061 


1.5 


3. 


105. 


GRB991105 


99 


6. 


57 


69 


0.07 


2.3 


6.8 


4.8 


GRB991213 


75 


5. 


55 


10 


0.05 


2.1 


9.3 


9. 


GRB000107 


115 


4.3 


90 


10 


0.05 


1.5 


3.3 


10. 


Mean 


93 


10 


60 


50 


0.064 


2.26 


9.4 


48 


Deviation 


13.5 


7.6 


9.6 


19.4 


0.014 


0.57 


5.0 


52 



